Search for: All records

This dataset originates from a new CESM2 CAM6 perturbed parameter ensemble (PPE) designed to explore climate and hydroclimate dynamics under a wide range of sea surface temperature (SST) conditions. The SST varies from 4 degrees Celsius colder to 16 degrees Celsius warmer than preindustrial levels, encompassing a broad spectrum of mean temperatures spanning the past 65 million years. This dataset offers valuable insights into climate and hydroclimate responses, as well as weather and climate extremes under diverse conditions.The dataset includes results from nine PPE simulations with different SST scenarios: preindustrial (PREI), 4K cooler (M04K), and 4K, 8K, 12K, and 16K warmer (P04K to P16K). For SSTs exceeding 8K warming, sea ice was removed to improve numerical stability. Each PPE set consists of 250 ensemble members, with 45 parameters related to microphysics, convection, turbulence, and aerosols perturbed using Latin Hypercube Sampling. An additional simulation with default parameter settings brings the total to 251 simulations, each running for five years using CAM6.3 (https://github.com/ESCOMP/CAM/tree/cam6_3_026; with additional paleo modifications).Post-processing converted the data into compressed NetCDF-4 format. All 251 runs were concatenated using ncecat to minimize the number of files. For example, the following file contains monthly surface temperature data from the preindustrial PPE: f.c6.F1850.f19_f19.paleo_ppe.sst_prei.ens251/atm/proc/tseries/month_1/f.c6.F1850.f19_f19.paleo_ppe.sst_prei.ens251.cam.h0.TS.000101-000512.ncA detailed variable list [https://rda.ucar.edu/OS/web/datasets/d651038/docs/detailed_vars.txt] can be found in the Documentation Tab.Parameter values are provided in the PPE Parameter File. More details can be found in the paper: Zhu et al. (2025). Investigating the State Dependence of Cloud Feedback Using a Suite of Perturbed Parameter Ensembles, Journal of Climate. 

Simulating the warmth and equability of past hothouse climates has been a challenge since the inception of paleoclimate modeling. The newest generation of Earth system models (ESMs) has shown substantial improvements in the ability to simulate the early Eocene global mean surface temperature (GMST) and equator-to-pole gradient. Results using the Community Earth System Model suggest that parameterizations of atmospheric radiation, convection, and clouds largely determine the Eocene GMST and are responsible for improvements in the new ESMs, but they have less direct influence on the equator-to-pole temperature gradient. ESMs still have difficulty simulating some regional and seasonal temperatures, although improved data reconstructions of chronology, spatial coverage, and seasonal resolution are needed for more robust model assessment. Looking forward, key processes including radiation and clouds need to be benchmarked and improved using more accurate models of limited domain/physics. Earth system processes need to be better explored, leveraging the increasing ESM resolution and complexity. 

Abstract. Climate variability in the last millennium (past 1000 years) is dominated by the effects of large-magnitude volcanic eruptions; however, a long-standing mismatch exists between model-simulated and tree-ring-derived surface cooling. Accounting for the self-limiting effects of large sulfur dioxide (SO2) injections and the limitations in tree-ring records, such as lagged responses due to biological memory, reconciles some of the discrepancy, but uncertainties remain, particularly for the largest tropical eruptions. The representation of volcanic forcing in the latest generation of climate models has improved significantly, but most models prescribe the aerosol optical properties rather than using SO2 emissions directly and including interactions between the aerosol, chemistry, and dynamics. Here, we use the UK Earth System Model (UKESM) to simulate the climate of the last millennium (1250–1850 CE) using volcanic SO2 emissions. Averaged across all large-magnitude eruptions, we find similar Northern Hemisphere (NH) summer cooling compared with other last-millennium climate simulations from the Paleoclimate Modelling Intercomparison Project Phase 4 (PMIP4), run with both SO2 emissions and prescribed forcing, and a continued overestimation of surface cooling compared with tree-ring reconstructions. However, for the largest-magnitude tropical eruptions in 1257 (Mt. Samalas) and 1815 (Mt. Tambora), some models, including UKESM1, suggest a smaller NH summer cooling that is in better agreement with tree-ring records. In UKESM1, we find that the simulated volcanic forcing differs considerably from the PMIP4 dataset used in models without interactive aerosol schemes, with marked differences in the hemispheric spread of the aerosol, resulting in lower forcing in the NH when SO2 emissions are used. Our results suggest that, for the largest tropical eruptions, the spatial distribution of aerosol can account for some of the discrepancies between model-simulated and tree-ring-derived cooling. Further work should therefore focus on better resolving the spatial distribution of aerosol forcing for past eruptions. 

Abstract. We examine results from two transient modeling experiments that simulate the Last Interglacial period (LIG) using the state-of-the-art Community Earth System Model (CESM2), with a focus on climate and ocean changes relevant to the possible collapse of the Antarctic ice sheet. The experiments simulate the early millennia of the LIG warm period using orbital forcing, greenhouse gas concentrations, and vegetation appropriate for 127 ka. In the first case (127ka), no other changes are made; in the second case (127kaFW), we include a 0.2 Sv freshwater forcing in the North Atlantic. Both are compared with a pre-industrial control simulation (piControl). In the 127ka simulation, the global average temperature is only marginally warmer (0.004 °C) than in the piControl. When freshwater forcing is added (127kaFW), there is surface cooling in the Northern Hemisphere (NH) and warming in the Southern Hemisphere (SH), consistent with the bipolar seesaw effect. Near the Antarctic ice sheet, the 127ka simulation generates notable ocean warming (up to 0.4 °C) at depths below 200 m compared to the piControl. In contrast, the addition of freshwater in the North Atlantic in the 127kaFW run results in a multi-century subsurface ocean cooling that rebounds slowly over multiple millennia near the Antarctic ice sheet. These results have implications for the thermal forcing (and thereby mass balance) of the Antarctic ice sheet. We explore the physical processes that lead to this result and discuss implications for climate forcing of Antarctic ice sheet mass loss during the LIG. 

Abstract. During the mid-Pliocene warm period (mPWP; 3.264–3.025 Ma), atmospheric CO2 concentrations were approximately 400 ppm, and the Antarctic Ice Sheet was substantially reduced compared to today. Antarctica is surrounded by the Southern Ocean, which plays a crucial role in the global oceanic circulation and climate regulation. Using results from the Pliocene Model Intercomparison Project (PlioMIP2), we investigate Southern Ocean conditions during the mPWP with respect to the pre-industrial period. We find that the mean sea surface temperature (SST) warming in the Southern Ocean is 2.8 °C, while global mean SST warming is 2.4 °C. The enhanced warming is strongly tied to a dramatic decrease in sea ice cover over the mPWP Southern Ocean. We also see a freshening of the ocean (sub)surface, driven by an increase in precipitation over the Southern Ocean and Antarctica. The warmer and fresher surface leads to a highly stratified Southern Ocean that can be related to weakening of the deep abyssal overturning circulation. Sensitivity simulations show that the decrease in sea ice cover and enhanced warming is largely a consequence of the reduction in the Antarctic Ice Sheet. In addition, the mPWP geographic boundary conditions are responsible for approximately half of the increase in mPWP SST warming, sea ice loss, precipitation, and stratification increase over the Southern Ocean. From these results, we conclude that a strongly reduced Antarctic Ice Sheet during the mPWP has a substantial influence on the state of the Southern Ocean and exacerbates the changes that are induced by a higher CO2 concentration alone. This is relevant for the long-term future of the Southern Ocean, as we expect melting of the western Antarctic Ice Sheet in the future, an effect that is not currently taken into account in future projections by Coupled Model Intercomparison Project (CMIP) ensembles. 

Variability in hydroclimate impacts natural and human systems worldwide. In particular, both decadal variability and extreme precipitation events have substantial effects and are anticipated to be strongly influenced by climate change. From a practical perspective, these impacts will be felt relative to the continuously evolving background climate. Removing the underlying forced trend is therefore necessary to assess the relative impacts, but to date, the small size of most climate model ensembles has made it difficult to do this. Here we use an archive of large ensembles run under a high-emissions scenario to determine how decadal “megadrought” and “megapluvial” events—and shorter-term precipitation extremes—will vary relative to that changing baseline. When the trend is retained, mean state changes dominate: In fact, soil moisture changes are so large in some regions that conditions that would be considered a megadrought or pluvial event today are projected to become average. Time-of-emergence calculations suggest that in some regions including Europe and western North America, this shift may have already taken place and could be imminent elsewhere: Emergence of drought/pluvial conditions occurs over 61% of the global land surface (excluding Antarctica) by 2080. Relative to the changing baseline, megadrought/megapluvial risk either will not change or is slightly reduced. However, the increased frequency and intensity of both extreme wet and dry precipitation events will likely present adaptation challenges beyond anything currently experienced. In many regions, resilience against future hazards will require adapting to an ever-changing “normal,” characterized by unprecedented aridification/wetting punctuated by more severe extremes. 

North Atlantic cooling during Heinrich Stadial 1 triggered an east-west precipitation dipole over the tropical Indian Ocean. 

Accurate understanding of permafrost dynamics is critical for evaluating and mitigating impacts that may arise as permafrost degrades in the future; however, existing projections have large uncertainties. Studies of how permafrost responded historically during Earth’s past warm periods are helpful in exploring potential future permafrost behavior and to evaluate the uncertainty of future permafrost change projections. Here, we combine a surface frost index model with outputs from the second phase of the Pliocene Model Intercomparison Project to simulate the near‐surface (~3 to 4 m depth) permafrost state in the Northern Hemisphere during the mid-Pliocene warm period (mPWP, ~3.264 to 3.025 Ma). This period shares similarities with the projected future climate. Constrained by proxy-based surface air temperature records, our simulations demonstrate that near‐surface permafrost was highly spatially restricted during the mPWP and was 93 ± 3% smaller than the preindustrial extent. Near‐surface permafrost was present only in the eastern Siberian uplands, Canadian high Arctic Archipelago, and northernmost Greenland. The simulations are similar to near‐surface permafrost changes projected for the end of this century under the SSP5-8.5 scenario and provide a perspective on the potential permafrost behavior that may be expected in a warmer world. 

Abrupt climate changes during the last deglaciation have been well preserved in proxy records across the globe. However, one long-standing puzzle is the apparent absence of the onset of the Heinrich Stadial 1 (HS1) cold event around 18 ka in Greenland ice core oxygen isotope δ 18 O records, inconsistent with other proxies. Here, combining proxy records with an isotope-enabled transient deglacial simulation, we propose that a substantial HS1 cooling onset did indeed occur over the Arctic in winter. However, this cooling signal in the depleted oxygen isotopic composition is completely compensated by the enrichment because of the loss of winter precipitation in response to sea ice expansion associated with AMOC slowdown during extreme glacial climate. In contrast, the Arctic summer warmed during HS1 and YD because of increased insolation and greenhouse gases, consistent with snowline reconstructions. Our work suggests that Greenland δ 18 O may substantially underestimate temperature variability during cold glacial conditions.